Effect of flow on the quasiparticle damping rate in hot QCD plasma 
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Abstract 

We derive an expression for the quark damping rate in hot quark gluon plasma in presence of flow. Unlike earlier 
treatments here all the bath particles are out of equilibrium due to the existence of non-zero velocity gradient. The 
quasi particle damping rate is found to be finite once the hard thermal loop (HTL) corrected gluon propagator is used. 



In recent years several authors have calculated the cmasrparticle damping rate (T) both in quantum chromodyamic 
(QCD) and quantum electrodynamic (QED) plasmas yj, y, l3i |4j, |5l |a |7[] . The calculations are important to understand 
various dynamical properties of such plasma associated with the quasiparticle excitations. Such investigations are 
important as the inverse of the damping rate gives the relaxation time used in the collision term of the Boltzmann 
equation. Therefore, once known, it can be used to estimate the thermalization time of a preequilibrium parton gas in 
relativistic heavy ion collision [8]. The expression of damping rate can also be used to derive the drag and the diffusion 
coefficients of the quasiparticle by folding the expression with relevant energy or momentum transfer iliolllillil 

[□HUB. 

For the relativistic hot plasma, it has been observed that quasi-particle life time or the damping rate calculations 
are plagued with infrared divergences. In non-relativistic plasma, where one considers only the coulomb or electric 
interaction such divergences are removed by the Debye screening effects. The problem becomes non-trivial in dealing 
with relativistic plasma where one has to worry about the electric and transverse interactions both. This additional 
complication actually arises due to the absence of static magnetic screening. In case of QED, in ((j |7[] the authors 
have shown that the electric contribution to the quasi-particle damping rate with plasma screening effects T is finite 
and is of the order r/ ons ~ g 2 T whereas the transverse part remains divergent. This is because the later is only 
dynamically screened. To obtain a finite result authors had to develop another resummation scheme to obtain finite 
result la, |7[]. In case of QCD one also encounters similar problem in hot plasma which however can be removed by 
assuming a 'magnetic mass' for the gluons. It would be worthwhile to note that in degenerate plasma one obtains finite 
results with magnetic interactions without further resummation unlike its high temperature counterpart. For detailed 
discussions about these issues we refer the readers to iHIt! Il6[ [TtIi . 

The purpose of the present work is to calculate the quark damping rate in QCD plasma with flow. Naturally this is 
a major departure from the above cited calculations where always the bath particles are assumed to be in equilibrium. 
This is true only when there exists no velocity or temperature gradient in the plasma and there is no external force. The 
scenario that we currently consider is different. Here, we assume that there exists a non-zero velocity gradient i.e flow 
in the plasma. Under such condition, one can no longer assume that the bath particles follow the ordinary equilibrium 
Bose-Einstein or Fermi-Dirac distributions for the interacting bosons or fermions respectively. The situation is very 
similar to what one encounters for the calculation of transport quantities like coefficient of viscosity or conductivity 
etc. We shall see in this letter that the transferred momentum dependence of the damping rate in presence of flow is 
very different leading to results which are finite once the plasma effects are included via the dressed propagator. 

To expose this further, we first recall that, the calculation of decay width in plasma is directly related to the collision 
integral of the Boltzmann equation given by the following equation: 



d \ 

^ + v p .V r + F.Vpj«p = -C[n p ], (1) 

here, p is the momentum of the quasiparticle, F is the external force. The right hand side of the equation is the collision 
integral. For the relativistic plasma, one assumes, that all the particles are moving at the speed of light, which means 
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that the velocity is a unit vector, v p = p = p/|p|. For two body interaction, which we assume here, the collision integral 
can be written as, 



r d 3 k d 3 p d 3 k' 

Clnp] 'J (2n)H2nfW) 



Xn p n k (l + n p ){\ ± n' k ) - n' p n' k (l + n p ){\ ± n k )(2n) 4 6 4 (p +k-p - k')\M\ 2 , (2) 



«, are the Fermi and Bose quasiparticle distribution functions for quarks and gluons and the ± includes stimulated 
emission or Pauli blocking. \M\ 2 is the relevant squared matrix element for the 2 — > 2 processes summed over final 
states and averaged over initial states. 

In absence of the external force and gradients of temperature, velocity or density Eq.(Q} takes very simple form. If 
n p is considered to be slightly out of equilibrium i.e n p - n p + 6n p , in relaxation time approximation the solution of 
the Eq.([TJ can be written as n p (t) - n p (Q)er r i> t , where T p is the damping rate of the quark moving with the momentum 
p related to the inverse life time t" 1 . Mathematically, 

C[n p ] = 5n p lT p , (3) 

where, t p is the relaxation time of the collision in this case and 1/t p = T p . The collision term C[n p ] consists of two 
parts. One involving equilibrium distribution function and the other non-requilibrium 5n p s, 

C[n p ]=C eil [n P ]+C" "- e nn p ]. (4) 

The C eq [n p ] vanishes due to the energy conserving delta function i.e, 

[n p n k (l -n p >){\ - n k >) - n p >n k >(\ - n p )(l - n k )]6(p + k - p -k'). (5) 

^fTnglpar 

Mm, 



From the C non e 5l^el£ art we 8 et tne following expression for the damping-rate when only one particle is slightly away 
of equilibrium | 



r<» = f [n k (l - n k >)(l -n p -)- n p >n k >(l - n k )\ (2n) 4 6 4 (p + k-p- k')\M\ 2 . (6) 

Jk,p',k' 

We use J k , instead of d 3 kd 3 p'd 3 k' /(2n) 9 throughout the letter. In the calculation of the above damping rate only the 
test particle is assumed to be slightly away from equilibrium while the bath particles are not. In this case we do not 
need to know about the explicit form of the non-equilibrium distribution function as it gets canceled from both sides 
in the Eq.(f3]l. On the other hand the general expression for the damping rate when four particles are out of equilibrium 

is, 

r(p) =11 ^ dn p (n k (l ±«/)(l ±n k >) +n p >n k >{\ ±n k )j 

+6n k (n p (\ ±n p >){\ ±n k >) + n p >n k >{\ ± n p )j - 6n p > \n k >{l +n p )(\ ±n k ) + n p n k {\ ±n k ')j 

-dn k - (n p -(l + n p ){\ ± n k ) + n p n k {\ ± n p >)) (2^)V(p + k-p - k')\M\ 2) j/6n p . (7) 

In the above equation it is evident that if we put 6n k = 6n p ' = 6n k ' = 0, then we get back the Eq.©. To obtain 
the expression of the Eq.((7]i we have to know the explicit form of (5n, (i = p, k, p',k ) unlike the above situation. The 
non-equilibrium distribution function takes different form depending upon the problem considered. For our purpose 
we take the distribution function as 



^l^C^^yX^ (8) 



s dpi 

where, O, A V = pi, x pi,yf(p/T) and C is some constant which has been estimated in 5Hl . In the presence of a small shear 
flow u(y) in the x direction X_ n , = du x /dy, f(p/T) is some rotationally invariant function depending only on the energy 
of the excitation which has to be determined from the variational calculation of the Boltzmann equation 1 1_8 , Ijl 2011 . 
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But a much simpler and standard way is to take the trial function in the viscous process as, f(p/T) = (p/T) 2 lfl9ll . 
Substituting the non-equilibrium distribution function Eq.© in the collision integral we obtain, 

Co» =H^ f n p n k (\±n p ){\+n k )\M\ 2 (2n) A 5\p + k-p -k) 

sl Jk,p,p' 

[(D, v (p) + <£, v (k) - 0,,(p') - 0, v (k')]. (9) 
To proceed further we introduce the following expression here, 

[GVGsJ =(-)^ f « / ,« t (l±«; ; )(l±n;)|M| 2 (27r)V( j p + ^-/ 3 '-^') 

[O^(p) + <D A ,,(k) - (D, v (p') - <D, v (k')] 2 . (10) 

We also define, 



3, 

sT 



[Q^Sxyl = i-)^n p {\ ± n p )f 2 (p/T), (11) 



where, S „ - 5n p . With the help of above two equations one can easily check that expression for the damping rate 
takes the following form, 

[t> AT ,CoJ 

r( " , = lOT <12) 

In the collision integral we use the spatial delta function to perform the k integration and to shift the p integration 
into an integration over q, where, q = p - p. Hence, 

r 0?) - FT I q 2 dqk 2 dk \ dcos6dcos6 k „ I dcpn k {\ + n ' k )\M\ 2 6(p + k - p - k ) 

(2n) p Jo J -i Jo 

[<D^(p) + O x> ,(k) - <£, v (p') - 0,,(k')] 2 . (13) 
It is convenient to introduce a dummy integration variable oj, 

da>6(to + p - p')S(u> + k- k). (14) 

oo 

Evaluating q = p - p in terms of p, q and co&9 pq and defining t — lo 2 - q 2 we find, 

P_ 

pq~\ m q 2pq 

k' I 0) t \ 

6{o + k-k) =—6\cos6 kq + — . (15) 

kq \ q Lkq j 

Now, to proceed further square bracketed term containing <b xy in the Eq.(fT3l has to be evaluated. Considering the 
small angle scattering limit which is a very reasonable assumption in the present context ll20ll . one can easily obtain 
the equation mentioned below. We need not worry about the cross terms since they do not exist after carrying out the 
integrations over da> and d(f>. So, the terms which actually contribute to the damping rate are the following, 

[<MP) - %y(p')] 2 = OJ 2 [f(p)'] 2 + 3^-^[f(p)] 2 , 

P 2 

[®^(k) - 0), y (k')] 2 = aj 2 [f(k)'] 2 + 3^-^[f(k)] 2 . (16) 

P 2 

Upto now, we have not made any reference to the interaction which is contained in the matrix amplitude. For QCD 
plasma, to calculate quark damping rate we need to consider qq — > qq as well as qg — > qg processes. For the former 
case the squared matrix element is given by, 

l^l 2 - y U ^A- (17) 



5(u> + p-p) — — 6\cos6 pq — 



The above matrix amplitude squared gives divergent contribution to the damping rate due to the presence of sin- 
gularity in r 2 . Hence, one needs to incorporate the screening effect. In case of small momentum transfer, q « p, 
k ~ T, from the energy conservation delta function we obtain a> = p — p - v p .q = -Vk.q. The velocity projections 

transverse to q have lengths |v p .t| = |Vk,Tl = V(l _ x2 )> where, x = co/q. Consequently, [Vp.xl-IVk.Tl = (1 - x 2 )cos(p, 
where, <p is the angle between the two velocity projections on q. For q « T we have, 



2pk( l - cos8 p k) 
2pk( \ - jc)(1 - cos^ 



(18) 



Hence, the expression in Eq.(IT7b takes the following form, 



\M, 



1 (1 



^ 2 ) 2 cos 2 ^ 



q 4 (q 2 - o?f 



(19) 



From the above expression it is evident that the matrix amplitude squared is infrared divergent. To circumvent this 
problem one adopts Braaten and Pisarski's method 1I22L In this formalism, the interaction is split into hard and soft 
momentum transfer domain by introducing an intermediate momentum scale q* 12311 . For the hard sector one then uses 
the bare propagator, for the soft sector, on the other hand, one loop corrected gluon propagator is used. To get the final 
result both contributions have to be added up and on addition the arbitrary momentum scale q* gets canceled. 

The dressed gluon propagator involves the longitudinal and transverse polarization functions which are given by: 



U L (q, at) 
U T (q, u>) 



x + 1 
x-1 

2- 



x x(l- x ) / x + 1 

— + In 

2 4 U-l 



(20) 



where, mo is the Debye mass. Further simplifications can be made if small angle approximation is made which gives 
tli - tn% and i{nlA)xm 2 D . 

With the dressed propagator for the soft sector the matrix amplitude squared takes the following form, 



\M q . 



(1 - X 2 ) 2 COS 



( q 2 + n L f (q 2 - u 2 + n T ) 2 



(21) 



It is to be noted that the cross term of the longitudinal and the transverse interaction vanishes after the d<p integration, 
hence, we omit the term here. 

Using the expressions of the Eq.([To*i> and the quark-quark scattering matrix amplitude square in the Eq.([T3~t for the 
electric interaction in the soft region we have, 



-A. 



5T S InT 
+ 



\%lnp 2 8I/7' 

5r 3 InT 5 
+ 



1tiT 5 \ C q q 3 dq 
8l7jJo (q 2 + m 2 D 



Slnp 2 81p 4 



log 



—I 



ml) 2 



(22) 



In case of the momentum region where q > q* we can ignore the term containing Debye mass in the denominator, 
hence, one finds for the damping rate, 



5T 3 



InT- 



Unp 2 %\p A 



log 



Qmax I 

q* I 



(23) 



On addition of the Eqs. (l22b and ( f23l the intermediate momentum scale q* gets canceled. Hence, we obtain a finite 
result in the longitudinal sector independent of the arbitrary momentum scale, 



I 



5T 3 



InT- 



\81tt/7 2 8I/7 4 
4 



log 



qmaxl 



(24) 



In case of the transverse interaction one can proceed in the same way like the longitudinal one and obtain the final 
result, 



,(P) 



-A- 



5T' 



\ \62np 2 162p 4 



log 



2q„ 



yfnniD* 15 



Hence, the final expression for the damping rate in presence of flow in the medium has the following form, 



\54jtp 2 5V 



l I ^ I , oi \ c imax\ I 

log — +21og - — 

' V^ 1 1 m D 1 JO 



(25) 



(26) 



From the kinematics the q„ mx here can be chosen as T From the final expression it is evident that at high 

temperature in presence of flow in the medium quasiparticle damping rate turns out to be finite and logarithmic in 
nature. Interesting to note here that the viscous relaxation time also exhibits similar logarithmic behaviour I5I [l8l [l9ll . 
This is not surprising as these quantities are related to each other upto an integral 02411 . 

To summarize in the present work we have derived an expression for the quark damping rate in QCD plasma in 
presence of flow. It has been shown that the quark damping rate in presence of a velocity gradient differs from the 
usual damping rate because of the presence of AV (p) terms. The latter, brings an extra co 2 term in the numerator which 
changes the infrared behaviour of F completely. Interestingly, here we find that the quasi-particle damping rate turns 
out to be finite without further resummation as mentioned in the introduction. It is important to note also that F in the 
present treatment is independent of the coefficient of viscosity 77 and the magnitude of the flow gradient, although the 
distribution functions contain both the terms. 
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